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Abstract 



In this paper we consider the relation between random sums and compositions of different processes. In particular, 
for independent Poisson processes Af„(t), Npit), t > 0, we show that N„{Np{t)) = ^ f i Xj, where the XjS are 
Poisson random variables. We present a series of similar cases, the most general of which is the one in which the 
outer process is Poisson and the inner one is a nonlinear fractional birth process. We highlight generalisations of 
vQ these results where the external process is infinitely divisible. A section of the paper concerns compositions of 

r — the form JV„(Tp, v e (0, 1], where zj, is the inverse of the fractional Poisson process, and we show how these 

QQ compositions can be represented as random sums. Furthermore we study compositions of the form 0(JV(t)), 

^SJ t > 0, which can be represented as random products. The last section is devoted to studying continued fractions 

, I of Cauchy random variables with a Poisson number of levels. We evaluate the exact distribution and derive the 

^^ scale parameter in terms of ratios of Fibonacci numbers. 

] Keywords: Fractional birth process. Bell pol5momials, Mittag-Leffler functions, Fibonacci numbers, continued 

. . fractions, Poisson random fields, golden ratio, Linnik distribution, discrete Mittag-Leffler distribution, logarithmic 

^ random variables, negative binomial distribution, Mellin transforms. 
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1 Introduction 

Publications in the field of probability have devoted considerable attention to compositions of different processes, 
as e.g. Brownian motions, fractional Brownian motions and telegraph processes. Also, more general cases as stable 
processes, combined in different ways have been investigated and the p.d.e. connections analysed. In the present 
paper we focus our attention on the composition of point processes, e.g. Poisson processes, fractional Poisson 
processes, and others. For independent homogeneous Poisson processes N^[t), Np[t), t > 0, we are able to show 
that JV„(JV^(t)) has a remarkable connection, with respect to distributional properties, with random sums, that is. 



we prove that 

Nfiit) 

JV(t) = JV„(JV/3(t))=2Xj.. (1.1) 

j=l 

The first part of the paper is devoted to the presentation of similar results involving more general processes as e.g. 
nonlinear fractional birth processes S^'^Ct), t > 0, v e (0, 1], so as to obtain 

9" it) 

where the XjS appearing in (1.1) and (1.2) are independent Poisson random variables of parameter A„. A more 
general result is obtained when the external process is replaced by a process possessing infinitely divisible 
distribution. In this case, we are able to show that 

Npit) 

where the random variables i^^s are the components of the infinitely divisible random variable 6(1), and Np[t), 
t > 0, is a homogeneous Poisson process. The representation (1.1) is remarkable in that it results in the explicit 
law of a random sum which is rarely possible in general. The distribution of (1.1) can be given as 

Pr{NjN^(t)) = k} = -^e-^^^2 ^^-^^ (1.4) 

~k\ 
where 



^e-^^'fi-^ ")2J,(A^te-^«), 



<B,(x) = e--2— p, (1.5) 



r! 

r=0 



are the so-called Bell polynomials, and A„, A^, are the parameters of, respectively, N^(t) and Np(t). The 
composition (1.1) produces a process with linearly increasing mean value and variance: 



|EJV„(JV^(t)) = A„A^t, 
\VarJV,(JV^(t)) = A„(A„ + l)A^t. 



(1.6) 



For the iterated Poisson process we find that the first-passage time 

T, = inf(s : JV„(JV^(s)) = k), k> 1, (1.7) 

has distribution 



Pr{r,eds} = dsApe-^«e-V^^e-^»^[a + l)''-;'=] -^, s < 0. 



(1.8) 



Furthermore, Pr{r^ < oo} < 1, for all fc > 1, so that there is a positive probability of never hitting the level k 
because the iterated Poisson process can take jumps of arbitrary integer-valued size. 



Finally, we note that the iterated Poisson process produces a Galton-Watson process (continuous in time) . 

In the case in which Np{t) is replaced by a non-homogeneous Poisson process 9T(t), with rate A(t), t > 0, we still 

have a representation of the composition JV„(91(t)) as a random sum, i.e. 



<yiit) 



N,mt))=Y,^,, 



(1.9) 



;=i 



where the XjS are independent Poisson random variables of parameter A„. Since the probability generating 
function of (1.9) reads 



Eu- 



N„(Sfi(t)) ^ J|jA(w)dw[e^»("-i'-l] 



we have also that 



fEJV„((t)) = Aj^'A(w)dw, 
lVar=A„(A, + l)rA(w)dw. 



(1.10) 



(1.11) 



The process (1.1) has non-decreasing sample paths with jumps of integer-valued size and thus differs and extends 
the classical homogeneous Poisson process. We note that, for A„ 7^ Aa, we have that 



Pr{jV„(JV^(t))#JV^(JV„(t))} = l, 



(1.12) 



as can be inferred from the structure of the probability generating function. The interchange of the non- 
homogeneous Poisson process 9T(t), with the homogeneous one produces a composed process 



^(iV„(t)), 
whose probabihty generating function can be written as 



^Kty 



^ rl 

r=0 



(1.13) 



(1.14) 



Even the mean value of (1.13) differs from (1.11) since 



E9a(JV,(t)) = e-^«^2 



r=l 



A(w)dw 



r! 






A(w)dw 



j-i 



Pr{JV„(t)>;}. 



(1.15) 



The third section deals with the inverse of the fractional Poisson process N^[t), t > 0. The process N^[t) can be 
viewed as a renewal process with Mittag-Leffler distributed intertimes and the following probability distribution 
(see Mainardi et al. [2004], Beghin and Orsingher [2009]) 



f^ ymj r(v; + 1 



(1.16) 



The inverse of the fractional Poisson process is defined as 

T^ = inf(t:iV^(t) = /c), fc>l, ve(0,l]. 



t >0, ve(0,l]. 



(1.17) 



with distribution 



Pr{Ti; e ds}/ds = A*^ > A!, (1.18) 

and moment generating function 

Ee-''^*=K- + 1 . (1.19) 

It should be noted that for v = 1 this coincides with the Erlang process. The composed process N^^tJ.) has 
probability generating function 



which suggests the following interesting representation: 



l + (l-iz)V" 



k 



(1.20) 



k 

JV„(tP^2?^, (1.21) 

j=i 

where the independent random variables ^j, I < j < k, are discrete Mittag-Leffler (see Pillai and Jayakumar 
[1995]). From (1.20), we can infer that N^(t^) has Linnik distribution and, for v = 1, this coincides with the 
negative binomial distribution having parameters k and Pia/i^a + ^p)- For l^he special case v = 1, we also have the 
following representation of JV„(Tp: 

JV„(t^)=Xi + ---+X^, (1.22) 

where JV is a Poisson random variable of parameter /x = logCCA^ + A^)/Ap)'^, and the XjS are independent and 
possess logarithmic distribution of parameter 17 = A„/(A„ + A^). Furthermore, the case iV„(</)^) where the inner 
process is the inverse 0^ of a fractional linear pure birth process YHt), is examined. 
In particular, we show that 

rf^^^i^ + i 



EuN„(y;(0) ^ fc, V -^ J |^| < j_ ^-^^23) 



r(Mlz!i): + i + fc 



The final section of this paper deals with random products of non negative i.i.d. random variables of the form 

N(_t) 



N„M = Y\Xj, t>0, (1.24) 



i=l 



where N[t) is a homogeneous Poisson process. We show that the Mellin transform of iV^(t) is 

E(JV„(t))''-' = e^^K""-i]. (1.25) 

Furthermore, we can evaluate the covariance function of the process N^[t) which reads 



Clearly, the random products can be viewed as compositions of the form S(JV(t)), where 

""=(n;..x,=.2,'.'«. III 

Finally, we consider continuous fractions with a random number of levels: 

1 

[^i>^2>--->^Nit)]=^i + ~ 1 > (1.28) 

X2 H 

where the XjS are independent standard Cauchy random variables and N[t) is a homogeneous Poisson process. 
We show that the conditional distribution of (1.28) is Cauchy in which the scale parameter equals F^^i/F^, where 
F„ are Fibonacci numbers. This permits us to give a stochastic representation of (1.28) in the form 

J^ w(t)+ l 
[^l',^2>--->^Nit)]= 2^ ^j,-N(t)' (1.29) 

where l^j,Ar(t) are Cauchy random variables with scale parameter equal to 1/Fjv(t). 

2 Composition of Poisson processes with different point processes 

In this section we examine the following compositions: 

1. N„[Np[t)), where JV„ and Np are independent Poisson processes; 

2. N„(9T(t)), where OT(t) is a non-homogeneous Poisson process with rate A(t), t > 0; 

3. iV„('3^^(t)), where '3^^(t) is a nonlinear fractional birth process with birth rates Aj, . . . , Aj,, . . . , and < v < 1; 

4. JV„(y^'(t)), where 7^' (t) is a linear fractional birth process. This is a special case of the preceding point. 

5. JV„(Jf(^)), where Jf(^) is a Poisson field. 

6. G(JV^(t)), where 9(1) is an infinitely divisible random variable and N[t) is an arbitrary point process. 
We will establish some distributional relations between these composed processes and random sums. 

2.1 Iterated Poisson process 

We start our analysis by considering the iterated Poisson process N^t) = JV„(JV^(t)), t > 0. The sample paths of 
JV(t) are non-decreasing with jumps of arbitrary integer- valued size. In the figures 1 and 2, we give the trajectories 
of JV„(t) and N^^t) separately and then the sample path of JV(t) obtained by their composition. 
We note that the iterated Poisson process N[t) jumps at the occurence of events of the inner process Np^t). Thus, 
if the rate of JVp(t) is large, then N[t) has rapidly increasing trajectories for large A„. However high values of A„ 
and low levels of An can produce contradictory results and thus compensate each other 
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Figure 1: A realisation of the external Poisson process N^{t), t > 0. 
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Figure 2: A path of the process JV(t) = JV„(JV^(t)), t > 0, together with its relation for a specific path of the internal 
process Np[t), t > 0. 



Theorem 2.1. The distribution ofN[t) = N^[Np(it)) reads 



K , ,^e"^"''rHA«t7 



where 



k\ ^ r! 

r=0 

= ^g-A,t(i-e-^«)25^ (Apte-^«) , fc > 0, t > 0, 



~ r'^x 



<B,(x) = e--^^:p, (2.2) 



is the kth order Bell polynomial [Boyadzhiev, 2009]). The probability generating function o/JV(t) has the form 

E^N(0 ^ gA,t(e^«(.'-i)-i)^ 1^1 < ^_ (2.3) 



Proof. 



Pr{iV(t) = fc} = Xi Pr{iV„(r) = fc}Pr{JV^(t) = r} (2.4) 

r=0 

00 

= 5'- fc! 



^_,„,(A^^_.,,(v):^ 



Furthermore 



e„...,,g;„.g;'-'-W.O>/-V(V)' ^^^, 



fci ri 

fc=0 r=0 



z 



^-A„r^LiA„r^A«t '- P ^ 



r=0 

and thus (2.3) emerges. D 

In Figure 3 we present the first four state probabihties as a function of time t for the iterated Poisson process in 
which A„ = A^ = 1. 

Theorem 2.2. The following equality in distribution holds: 

NSNpM)=Y,Xj, (2-6) 

;=i 

where the XjS are Ltd. Poisson random variables of parameter A„. 
Proof. The probability generating function of the random sum (2.6) is 

Eu^i=i ""' = Y, (Eiz^)' Pr{iV^(t) = fc} = 6-^^^+^^^="" = e^^^('='"'"""-i), (2.7) 

k=0 

and this coincides with (2.3). D 

7 
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Figure 3: The first four state probabilities (Iterated Poisson). The parameters are A„ = 1, A« = 1. 



Remark 2.1. For the iterated Poisson process we can write 

EiV(t) = EiV,(iV^(t)) = A„A^t, 



(2.8) 



YarN(it) = A„(l + A„)A^ t. 



(2.9) 



We can obtain (2.8) directly, by means of the probability generating function, and also by applying Wald's formula for 
random sums. 



00 oo 



k=0 r=0 



EN„(iV^(t)) = ^fc2e 

r=0 

oo 



XKr)\.,^MptY 



k\ 



r! 






e "'-^^e 

r=0 



r! ^(fc-1)! 



-Xf,t 



/.^a '- ~ ^r,^Rt. 



r\ 



Alternatively 



EN„(JV^(t)) = -Eiz*^^) 



d 
du^ 



A„A^t(e^«f"-i))e^^^C'='"'""'") 



■ A„A«t 



u=l 



and, by applying the Wald's formula 

EXEJV^(t) = EJV(t) = A„A^t. 
For the second moment, analogously we obtain 



00 oo 



EJV' 



^co=i:fc^i:e- 



AKrf.Sh^Y 



k=0 r=0 



k\ 



r\ 



(2.10) 



(2.11) 



(2.12) 



(2.13) 



s^ 



e-^^^> e-^' 



r! fH fc! 



Sm(^«'-)^ 



fc=o 



= e"^'^' 



r=0 ■ J:=0 



S^ 



e-^^'> e-^' 



,a^0'" 






k=0 



k=0 



k\ 






= e"''-^' 



A„A^t + A^A^t + A^A„t . 



Therefore 



VarJV(t) = EJV^Ct) - [EJV(t)] ^ = A„A^(1 + A„)t. 
This result can he confirmed by applying the Wald's formula for the variance. 

Yarffit) = YarN/^it) [EX]^ + YarXENp^t) = A„A^(1 + A„)t, 
because EX = YarX = A^. 
Remark 2.2. The state probabilities pj^it) = Pr{JV(t) = k} satisfy the difference- differential equations 



'' A!" 



dt 



Pk 



(0 = -^pPkM + A^e-^« Y, -^^Pk—^^^' ^-^ 



In order to prove (2.16) we write 



p,(t) = Y,Pr{N^[r) = k}Pr{Np[t) = r}, 



r=0 



SO that 



J 00 J 

^Pfc(t) = 2Pr{JV„(r) = k}-Pr{Npit) = r} 

oo oo 

= -Xp Y.P'-i^air) = k}Pr{Np[t) = r} + Xp ^PrjJV^Cr) = fc}Pr{JVp(t) = r - 1} 

r=0 r=0 

= _ A,p,(t) + A, ^e-^.^ V e-' y}^''\~' 
pfkK ^ P^i ^ (r-1)! 



r=0 

-A^i3,(t) + A^e-^^^-fe-^' 

' ' fc! -^m!(_fc — mjl 



fc! 



y — - — y 



g A„, J.,, 



(2.14) 



(2.15) 



(2.16) 



(2.17) 



(2.18) 



m=0 



k ^m 



m=0 

It is now easy to show that 

00 

G(iz,t) = 2iz''M0 (2.19) 



k=0 

satisfies the partial differential equation 



— G(u, t) = -A^G(u, t) + A^e^«("-i^G(iz, t) = A^G(iz, t) [e^-^"-^' - l) . (2.20) 



Remark 2.3. For the composition JV„(9T(t)) we have that 



00 00 



r:;i r^'^ ,i f A(w)dw 
X;"'MiV„(9T(t)) = k} = J^u^Y^^-^e-fo AWd«'i:^^_ Le-A„r (2.21) 



fc=0 J;=0 r=0 



^ rrA(w)dwi'' 

V e-^-'^ ^e^«™e- ^0 ^W^h 

/-J rl 



r' 

r=0 



^^[eAa("-i)_i]JjAWdlv_ 

We can aZso ascertain that 



m(_t) 

d 



iVam(t))=Xi^i (2.22) 



by using result (2.21). 



2.1.1 Hitting time for the iterated Poisson process 

Here we study the distribution of the random variable 

T^ = inf(s : N,[Npis)) = fc), fc > 1, (2.23) 

which represents the first-passage time of the iterated Poisson process at level fc. In the next theorem we state the 
main result. 

Theorem 2.3. 

-ifc 00 /--I y 

Pr{rfceds} = dsA^e-^«-^e-V2e-^«^"[0 + l)'^-/]^^, s > 0. (2.24) 

j=o ■' ■ 

Proof. In order to arrive at result (2.24), we first write 

k 

Pr{r, e ds} = 2Pr{JV„(JV^(s)) = fc - h,JV„(JV^(s + ds)) = fc} (2.25) 



h=l 



10 



Y,Pr{N^iNp[s)) = k- h,JV„(iVp(s) + dJV^(s)) = k}. 



h=l 



Clearly, dJV^(s) either takes the value with probability 1 — A^ds (and, in this case, all events appearing in (2.25) 
are mutually exclusive) or the value 1. Therefore 

it 
PrlTfc e ds} = A^ds^PrfiV^CiV^Cs)) = k- h,N^[Np[s) + 1) = k}. (2.26) 

h=l 

In the interval (iV^(s), JVp(s) + 1), the external process can take all possible values 0<fc — h<fc — 1. The value k 
must be excluded because if at time s, iV„(JVp(s)) = k, the first attainment of value k cannot be recorded during the 
interval [s,s + ds]. Furthermore 

Pr{JV„(JV^(s)) = k- h, JV„(JV^(s) + 1) = fc|JV^(s) = ;•} (2.27) 

= Pr{JV„(;) = fc-h,JV„0 + l) = fc} 
= Pr{JV„a) = k- h}Pr{JV„a + 1) - JV,0) = h} 

(fc-h)! h\' 

By inserting (2.27) into (2.25), we arrive at 



{k-h)\ h\ ;! 



Pr{r,.ds} = dsA,e-V22e-»^-^^;^e-»9^4^ (2.28) 

j=0 h=l 



'^- ;=0 }■ 



We note that Vx{T^. < oo} < 1, for all fc > 1. From (2.24), we have that 



■yk 00 

Pr{r, < 00} = e-^«-^2e-^«^' [(; + l)*^ - ;''] (2.29) 



jk 00 /-j+l 



k 00 (-J+i 

J- j=0 J J 

^£ e-^»(^«V-idx 



(fc-1)! 

k oo c j+l 

(fc" 



< ^ I e-^-^x'^-Mx = 1. 

Cfc-l)!Jo 



D 



Remark 2.4. The previous result shows that there is a positive probability of never reaching level kfor the iterated 
Poisson process. For some cases this probability can be evaluated explicitly. 

Pr{T-i^ < 00} = A„e"^«2_is~^°'' = " _^ < 1- (2.30) 

;=0 1 ~ ^ 

11 



This is because 

< 1 - e-^» - A„e-^» = 1 -Pr{JV„(l) = 0} -Pr{JV„(l) = 1}. 
From (2.25) we can evaluate the distribution of Tj as follows. 

PrlTi e ds} = ds A^e"^- A^e-^P'^i-'="'°^, s > 0. 
By similar calculations we have also that 



2„-A„ 



Pr{T2 & ds} = ds ?i« 



A^e 



.g-Ap.(i-e-^«) |^i + 2(A^s)e-^«] . 



Furthermore 



Pr{T2 < 00} = [PrjTi < oo}] + —PrlT^ < oo}. 



Finally, an alternative form of (2.25) reads 



00 rj+i 
Pr{r;,eds} = A^ds^ e 

;=o h 



-KU+l-x)p^^ja^^s}, 



(2.31) 



(2.32) 



(2.33) 



(2.34) 



(2.35) 



where T^ = inf{s : N^[s) = k} and Pr{r" e ds} is the Erlang distribution for the external Poisson process. 

2.2 Subordination of a homogeneous Poisson process to a fractional pure birth process 

Let '3^^ it), t > 0, < V < 1, be a nonlinear fractional pure birth process with rates Xj > 0, j e N (representing an 
extension of the classical nonlinear pure birth process). It has been shown in Orsingher and Polito [2010] that 



Pr{3^^(t) = fc|3'^(0) = l}: 



1 lj=l ^j Z-im=l Til' 



E,,A-Kt") 



,.a,-A„) 



, k>l, 
fc = l, 



(2.36) 



where E^ i(x) is the Mittag-Leffler function. For the process N^iW^^tJ) we have the following result. 
Theorem 2.4. For the composition N^i'WXt)) we have that 

where the Xp j > 1, are i.i.d. Poisson random variables with rate A„. 

Proof. In order to prove (2.37), we evaluate the probability generating function of both members. 



(2.37) 



Eu^«(^''W^ = ^iz'' 



k=0 



— -E,,i-x,tn+ye-^'''-^r\x,y—, 

r=2 ; = 1 m=l 1 l; = l,;^m^^' ^mJ 



k\ 



Ev,li-^mi^') 



(2.38) 



r-l 



LXJ I — !_ I -J 

j=l m=l 1 l( = l,;^m(-^; ~ ^m) 



r=2 



12 



The probability generating function of the right-hand side of (2.37) reads 

Eu^r=i"'^i =e^«^"-i)pr{3^^(t) = l} + 2e^»f"-i^''Pr{3^^(t) = r}, (2.39) 

r=2 

as Eiz'^ = e^''~"~^\ because X is a Poisson random variable of parameter A„. By comparing (2.38) with (2.39), in 
view of (2.36), we have the claimed result. D 

Remark 2.5. For the fractional linear birth process Y^(t), t > 0, the distribution (2.36) specialises and takes the 
form 

PriY^M = fc|7-(0) = 1} = 1] f'^I J] [-iy-'E,^,[-Xpjtn. (2.40) 

For V = 1, the distribution (2.40) reduces to the geometric distribution of the classical Yule-Furry process. The 
probability generating function of the composed process JV„(7^(t)) becomes 

E[.-«(''^W)]=2.''2e-^«^^2(':_')(-l)^-%,(-A,;f') (2.41) 

k=0 r=l j=l ^^ J 

r=l j=l ^-^ ^ 

= 2(-l)^-Xa(-A,;t^)2 ( ': _ J )e-^«(i-"> 

oo oo ^ • -N 

= 2(-l)^-i£,,i(-Ap ;t^')e-^«(^-"^^' Xi (-l)'-e-^«(i-"> 

j=l r=0 ^ *" >^ 









j=i 



1 _ e^a("-i) 



2.3 Subordination with a Poisson field 

Let B a Borel set and let us indicate with A(B) its Lebesgue measure. The aim of this section is to study the process 

Nai^B)), (2.42) 

where N^{t), t > 0, is a homogeneous Poisson process with rate A^ > 0, and ^(B), is a homogeneous Poisson field 
with rate A > 0. Analogously to the iterated Poisson process, the representation 

iV„(^(B))=Xi + ---+X<^(B), (2.43) 
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where the random variables Xj are i.i.d. Poisson distributed of parameter A„ > 0, is still valid. Therefore, the state 
probabilities can be written as 

Pr{iV„(J?(B)) = fc} = -^e-^'^(^»-'="'°^'Bfc(AA(Be-^«)), fc > 0, (2.44) 

and the probability generating function reads 

Giu,B) = e^A(B)(.^»C'-')-i)^ 1^1 < ^_ (2.45) 

The emptiness probability is given by: 

Pr{iV„(J2^(B)) = 0} = e-^'\(^)(i-'="'""). (2.46) 

Let now B; be the disc with radius Z, centred in the origin. The first-contact distribution H(l), I e R"*", (with the 
first point) is 

H(Z) = l-Pr{JV„(J2^(B,)) = l} = l-e-^''''(i-'="'°\ ZeK+. (2.47) 

The probability density is in turn 

Zi(Z) = 2A7rZ(l-e-^«)e"^''''^^"'"''"\ Z e R+, (2.48) 

that is, a Rayleigh distribution (see, for the classical non subordinated case, Stoyan and Stoyan [1994], page 213). 

3 Compositions of Poisson processes with first-passage time of different 
point processes 

In this section, we consider a fractional Poisson process N^[t) (see Beghin and Orsingher [2009] for information 
on this process) whose first-passage time 

T^ = inf(t:JV^'(t) = fc), fc>l, (3.1) 

is composed either with a homogeneous Poisson process N„(t) (A„ > is its rate) or with a Yule-Furry process 

The distribution of (3.1) has the following density 

Pr(T^ e dt)/dt = -Pr{JV^(t) > fc} = A^t^"^-i£^^_^,(-Apt^), t > 0, A^ > 0, (3.2) 

[Beghin and Orsingher, 2010, formula (1.6)], where 

^'^ ■^r(^r-|-7)r! 

The function E^Sz) is called generalised Mittag-Leffler function (see Mathai and Haubold [2008, page 91]) 
Theorem 3.1. The composed process N^[k) = JV„(tP has the distribution 

1 ^ r-fcA f?^A ^^' r(v(fc + j) + r) 
Pr{N^[k) = r} = - y T^ ' „ ^ ,, , r > 0, (3.4) 

rij^V; J {kJ r(v(fc + ;)) 
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and possesses probability generating function 



2 '^' 

i+(i-izr-^ 



Iu|<l. 



Proof. We start our proof with the following relation: 



/* 00 

Pr{jV^(fc) = r}= Pr{JV„(s) = r}Pr{T^eds}, fc > 1, r > 0. 
Jo 



(3.5) 



(3.6) 



Instead of writing the distribution of t^ in terms of generalised Mittag-Leffler functions, it is more convenient to 

work with the following expression which we derive from scratch. 

From Beghin and Orsingher [2009], we know that the distribution of the fractional Poisson process is 



Pr{JV-(t) = m} = VC-D^-™ L /.^-, 



(Apt^'rX] 



; + m 



1) 



^V ;• jr(v(m + ;)+l)' 



t >0, ve(0, 1]. 



By considering that PrJTT < s} = Pr{JV^(s) > fc}, we have that 



PrfTj; e ds}/ds = ;r XI P^i^^C^) = ^} 



ds 

d 

ds 



m=k 



m=k j=0 ^ J J "^ \- 

dOO 00 y' ■ I ,„>. 



v(m + ;•) + !) 



h=j+m d 

ds 



ds 






■H^-^A 



^vh 



s'" i-fh-l 

-C-1) 
dsf-^^ "'"' r(vh+l) lfc-1 

h=k y~ y \. 



oo 



h— 1\ vhs 



^vh-l 



h=k 



fc-ijr(vh + i) 



j=h-k 



Xi(-A^y+'^(-i) 



fc - 17 r(vh) 

■j + fc-lA s^^'^+J)-! 



j=0 



fc-1 jr(v(fc + ;)) 






^v()c+j)-l 



'r(v(fc + ;)) 



(3.7) 



(3.8) 
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In the fifth step of (3.8), we used the following formula which is interesting in itself 



i]n(-ir=(-i)^f^:; 



(3.9) 



We provide here a proof of (3.9) in the following way: 

A ^, ^ h ^T ^ Jc-1 



m=k ^ ^ m=0 ^ -' m=0 ^ ^ m=0 ^ -' 

h{h - 1) hQi - l){h - 2) hQi - l)(h - 2)(h - 3) 
= -1 + h 1- - 



+ ••• + (-1)'^ 



2 2-3 2-3-4 

h{h-l){h-2)...{h-k + 2) 

2-3---(fc-l) 



h h{h - 2) hQi - 2){h - 3) 
~ 2 2-3 2-3-4 



+ ••• + (-1) 
(h-l)(h-2) 



^h(h-2)...(h-/c + 2) 



2-3---(fc-l) 



-1 + 



h 
4 



(h-l)(h-2)(h-3) 

2-3 
[h - l)(h - 2)(h - 3) . . . (h - fc + 2) 

(fc-2)! 
(h - l)(h - 2)(h - 3) . . . (h - fc + 2) 



h h{h - 3) 
3 3-4 

,h(h-4)...(h-fc + 2) 
!__ + .. . + (_l)fc 



,h(h-3)...(/i-fc + 2) 

+ • • • + (-1)'^-^^ — !^ 

^ ^ 3-4---(fc-l) 



4-5---(fc-l) 



{-lf-' + {-ir 



fc-1 



(-1)' 



(fc-2)! 
(h-1)! 



(-1)' 



-1 + 



fc-1 



(fc-l)!(h-fc)! ~*^~-^^ l^fc-l 



By inserting (3.8) into (3.6) we obtain that 
Pr{JV^(fc) = r} : 



-A„s 



iKsJ 



%{:r^ 



Mk+i)-\ 



a 






r(v(fc + ;)) 
r(v(fc + ;) + r) 



-ds 



I r+v((c+j) 



1 

r! 



j=o ^ •' 



X" 



k+j 



This proves formula (3.4). The probability generating function 



r=0 Jo 



°° -A 



tA.s)'- 



r! 



^^g(i>'* 



g-Vg^^ 



;=0 ^ •' 



r(v(fc + ;)) 

r(v(fc + j) + r) 

r(y(fc + ;)) • 



^v{k+j)-l 

r(v(fc + ;)) 

v()c+j)-l 

-ds 



ds 



r(y(fc + ;)) 



(3.10) 



(3.11) 



(3.12) 
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-")V § 



Afi 



A!(l-izr 



1 + 



-]0 


f h 




j J 


h V 


-k 



K^i-uy 






D 



Remark 3.1. The waiting time t^', fc > 1, for the kth event of the fractional Poisson process can be viewed as the sum 
of independent waiting times t\ separating the events of the Poisson flow, i.e. 






;=i 



It is well-known (see e.g. Mainardi et al. [2004] or Beghin and Orsingher [2009]) 

Pr{Tl.&<is}=Pr{Tl&ds} 



v^'-^S 



j=0 



rCvO + 1)) 



d 



ds 



Ev,li-^li^^)- 



From (3.14) it follows that 



PriTl- e ds} = 1, 



and, by writing the distribution of t^ as convolution of the terms pertaining to t^' -, we have also that 

Pr{Tl&ds} = l. 



The Laplace transform of (3.8) is easily calculated and reads 



Ee-'^^t 



-fc 



A*" °° 
_^ V 

IjVk Z—k V , 






^v(t+j)-l 

'r(v(fc +;•))• 



An 



-ds 



and this clearly confirms the additive structure (3.13). 

Remark 3.2. Result (3.17) suggests the following representation for the composed process iV^ (fc); 

iv^(fc) = ?I + ••• + ?!, 



(3.13) 



(3.14) 



(3.15) 



(3.16) 



(3.17) 



(3.18) 
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where the E,^. are independent random variables which are called discrete Mittag-Lejfler random variables (see Pillai 
and Jayakumar [1995]) having parameters v and X^/Xp. The discrete Mittag-Lejfler reduces to the geometric random 
variable of parameter q = A„/(A„ + Xp) when v = 1. We now give the explicit distribution of the generalised geometric 
random variables E,^'. 



PrW 



■'^§(1) 



j+i 



r(vO- + i) + r) 
rCvO + 1)) 



(3.19) 



Kr 


Jo 


h 


f OO 


^^J 






|^r(vO- + i))V K 



dw 



w 



r\ 

For V = 1, we extract from (3.19) the geometric distribution. 

, Xb f w'' ^-i)^ ,,, 

Pr{^i = r} = -^ e"*^— e"^"dw : 



-J-w" Idw. 
A! 



WW 



:pq' 



(3.20) 



where p = A^/(Aq + A^). In order to check that the generalised geometric law (3.19) sums up to unity, we write 

Xr i , / A 

Pr{^^ = r} = 

r=0 



r=0 a Jo y. '^aj 






dw 



-^v.l - 



(3.21) 



A/; 

A! 



1. 



Furthermore, formula (3.5) shows that N^{t} possesses Linnik distribution and its form is explicitly given by (3.4). 
Finally, the distribution (3.4), for v = 1, becomes the negative binomial distribution having parameters k and 
A„/(Aq + A«). Indeed, from (3.4), we have that 



Pr{N\k) = r]=Y,i-^r 



m=0 

'k + r-1 
r 

k + r-1 
r 

k + r-1 
r 

k + r-1 
r 



X.V^'^ik + m + r-l)] 



S 



m=0 \ a 



r!m!(fc-l)! 

k+m 



(3.22) 



(-1)" 



k + m + r — 1 



m 



S 



ay m=0 

k 



1 + 



-(fc + r) 
m 



A„ 



-(fc+r) 



A„ 



. Aa + A^ y \^ A„ + A^ 
Remark 3.3. We now find the distribution of a slightly modified first-passage time 



(D 



'^t' 



(3.23) 
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where t"^ is defined in (3.1) and has distribution (3.8). 



Pr{fl^ds} = ds(-Y\iy( ^Vli ., .. 

' VtJ PjL^yj J P r(v(fc + ;)) 



v[k+j)-l 



(3.24) 






The Laplace transform of (3.24) becomes 



r(v(fc + ;•))■ 



e-'''Pr{f^' e ds} = 



1 + A,- 



For /c — > cx) we obtain the fine result 



lim e-''^Pr{f^eds} = e-^'^'''". 



Result (3.26) shows that the rescaled first-passage time (3.23) converges in distribution to a positively 
law of order v e (0, 1). 



(3.25) 

(3.26) 

f skewed stable 



We now consider the Yule-Furry process Y^{t), with a single progenitor, subordinated to the first-passage time t^. 
The distribution of Y^{t^) is given below and can be determined as follows. Bearing in mind the distribution (3.8), 
we have that 



f°° ^ r-k^ 5v(«:+;)-i 

Jo 7=0 



^0 h=l 



(ifA^M::<i) 



the distribution (3.8), 
(3.27) 



J 

■^ (A„h)^('^+j:i 



h=l 



Jyh- 

a 

r-l 



i:ll<-"'- 



1 + 



Afi 



A^h^' 



a;; 

Ar 



The probability generating function of the distribution (3.27) reads 

r=l h=l ^ ^ 



r >0. 



>,h-i 



OO 



h=\ 



A^ 



-1 — t- 

a;^ 

1 + h''-^ 

Aft 



rr-1 



(3.28) 
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S-'"-(r^)' 



i + ^V 



Iu|<l. 



We remark that the inversion of sums in (3.28) is vahd only for |iz| < 1. 

3.0.1 The classical case 1/ = 1 

For V = 1 we have special interesting results for iV„(T^) and 7„(t^). For the first process we have the following 
result. 

Theorem 3.2. The composed process iV„(T^) has the following representation: 



JV„(t1)=Xi + ---+Xn, 



where N is a Poisson random variable of parameter 



M = log 



Ac 



(3.29) 



(3.30) 



and the XjS are i.i.d. random variables with logarithmic distribution of parameter q = A„/(A„ + A^). 
Proof. The random variable JV„(t^) is a negative binomial W (see (3.22)) with distribution 

'k + r-r 



Pr{W = r} : 



p'q'- 



(3.31) 



In our case p = A^/(Aq + A^) and q = Pia/i^a + ^p)- It is well-known that it can be expanded as a random sum of 
the form 



N,[tI)=X^ + ---+X^, 



(3.32) 



where iV is a Poisson random variable of parameter /i = — fc logp and X is a logarithmic distribution of parameter 

q. n 



Remark 3.4. From (3.22) we can infer that 



and 



EiV(fc) : 



A„ + Afl 



Af 



-k, 



(3.33) 



YarN(k) = —k. 



(3.34) 



These results can be confirmed by applying Wald's formula to the random sum (3.29). 
Remark 3.5. For v = 1, the distribution (3.27) becomes 

-k 



Pr{Y,^Tl) = r}=Y,f[_]\-ir' 

h=l ^ ^ 



l+h 



Afl 



r > 1, fc>l. 



(3.35) 
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We are able to give a fine expression of (3.35) /or fc = 1. We have 



P.fl'„(.;) = .! = g(::;)c-I)'-;^ C3.36) 



grn<-«"'' ' 



h J" 'A, 
In light of the formula (see e.g. Kirschenhofer [1996]) 



n^o^ ■' ^ ""(l + ^+^)' 



k=0 

the probability (3.36) becomes 



g(J^-^^^^ = x(x+i)-(x+^r ^^-^^^ 



Pr{y,(T ) = r} = / ,, ^ ° . ^ =/Beta r,/ + l , r > 1. (3.38) 

A„ r(^j.ij.^> A„ \ A„ ' 



e+'+o 



"a r n^ + 1 + r 



We can easily check that (3.38) sums up to unity because 

V^Betafr,^ + l')=V^ f 'x^Hl - x)VA„dx (3.39) 

^a Jo 

Remark 3.6. The mean value and the variance o/7„(t^) can be obtained by means of the following calculations. 

EYJtI) = ^— e^«^s'^-ie-Vds = ^— e-^^f>-^-^'s''-^ds = ^— , (3.40) 

" ' (fc-i)!Jo (fc-i)!Jo yh-Kj 

if X^ > X^. Analogously we have that 

^"'■^"^^^^ = W^. I ''^'^' + eV).-ie-Vds = i^-^J + (^^-^ J , A, > 2A„. (3.41) 

3.1 Composition of Poisson processes with the inverse of an independent fractional 
linear birth process 

Let YHt), t > 0, be a fractional linear pure birth process with rate A^ > 0, studied in Orsingher and Polito [2010]. 
From Cahoy and Polito [2010], the distribution of 

(t>l=mt:Y^[t) = k), (3.42) 

is obtained and has the following probability density: 

Mr,^<it}/dt = Y,f]f'^~^\-iy-'Xplt--'E,^,[-Xpltn (3.43) 

m=i (=1 V ' -'■ y 
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k ^ /' 1 \ 

1=1 m=l^ J 



t >0, ve(0, 1]. 



The relation 2m=; (T-i ) ~ (/)' ^^^^ iii t^he second step of (3.43), can be proved as follows: 



m=l ^ 



2-3 



2-3---(fc-Z) 



a + 1) 
a + i)a+2) 



z za+2) za+2)...(fc-i) 
1+ - + ^ — ^ + --- + ^ — - — ^ — - 

2 2-3 2-3---(fc-Z) 



Z Z(Z+3)...(fc-l) 

^+3+-+ 3-4-(.-Z) 



(Z + lXZ+2XZ+3)...(fc-l) 



(fc-Z-1)! 



1 + 



Z 



fc-Z 



The distribution of ZV„((/)n therefore becomes 



Pr(JV„(</,n = r) : 



°° „-A, 



e-^^\K^Y 



• (=1 v-y n=o 



^pM r(v(n+l) + r) 
'a; J r(v(n + l)) 



The probability generating function of iV„(0n can be written in a neat form as 



r=0 Jo 



~ '■~ (A„sy '^ 



r! 



2 [^)(-iy-'A^Zs^-i£,,,(-A^Zs^)ds 



[ e-^«^(i-") ^ (^^ j (-l)'-iA^Zs^-i£,,,(-A^Z5')d5 



;=i 



S <-)" 



A«Z 



Jc-l 



[A„(l-iz)]VApZ ^ 
fc! 



^s(':'V-«'^ ' 



A« 



+ 1 + Z 






A'.(l-u)'' 



+ 2 



A-d-u)" 



+ fc 



r(^IM: + i + fc 



fc • Beta fc, — 



Vii-uY 



Xa 



+ 1 , |IZ|<1. 



(3.44) 



(3.45) 



(3.46) 
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4 Poisson random products and Poisson random continued fractions 

4.1 Multiplicative compound Poisson process 

In this section we consider a multiplicative compound Poisson process (denoted here 7i-compound Poisson process), 
defined as 



N(t) 

N^M = Y\Xj, t>0, 



(4.1) 



where the XjS are i.i.d. random variables and N(t), t > 0, is a homogeneous Poisson process with rate A > 0. We 
start by calculating the Mellin transform oiN^^t). 



E [jV,(t)] "-1 = Y, [eX^-i] ' ^e-^^ = e^^f^""-!). 



k=0 



The relation (4.2) can be rewritten as 



(4.2) 



(4.3) 



For the non-negative random variables XjS, the random sum ^ ■ = l'^'^'^ logXj can be reduced to a Poisson random 
product for the random variables XjS possessing Mellin transform at point 17 = i^ + 1, 



EX''"^ = EX'^. 
We give the explicit form of the covariance function in the next theorem. 
Theorem 4.1. For < s < t, the covariance of the random product N^[t\ t > 0, reads 



(4.4) 



sTEXiX-l) 
At(EX-l) fpAsE[J>f(X-l)] _ pAs(EX-l)] ^ pAt(EX-l) ( Xe^^dw. 

sE(X-l] 



Cov(iV,(t),iV,(s)) = e^"-^-" [e^^=[^(^-"] - e'-^'^-'^] = e 



Proof. 



E 



■jV(t) N[s) 

j=l r=l 



E 



'N[s) Nit) Nis) 

n^r n ^rfi^. 

j=l l=Nis)+l r=\ 



(4.5) 



(4.6) 



22 [EX2]'"[EX]"-"'Pr{JV(s) = m,JV(t) = n} 

00 00 

2 2 ['^^l "" [lE^]"'™ Pr{iV(s) = m}Pr{JV(t - s) = n - m} 



m=0 r=0 

As(EX2-l) A(t-s)(EX-l) 



-A(t-s) 



(A(t-s))'- 



Therefore 



;ov 



Nft) N(s) 

.j=l r=l 



As(EX2-l) A(t-s)(E)<'-l) _ At(EX-l) As(EX-l) 



(4.7) 
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At(EX-l) [ ^XsE[XliX-l)] _ ^AsCEX-l)" 



rgA5E[X(X-l)] _ gMEf-l)] _ 
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Remark 4.1. Formula (4.5) shows that the process N„[t) is positively correlated. 
As a consequence of the previous calculations we have that 

and 

YarN^^t) = e'^^^'^'-i) - e2^fO^-i5 = ^-^ti^-^'') h _ g-AtECx-D^l _ ^^^^ 

For X ~ iV(0, 1), the covariance function ofN^(t) takes the form 

CovJV^(t) = 2sinh [Amin(s, t)] 6"^"'"^''^) = 2sinh [ACov(JV(t),iV(s))] e-^cov(jv(t),N(.))_ ^4 ^^q) 

Remark 4.2. if the random variables Xp j > 1, are positively skewed stable with index v e (0, 1), we are able to give 

an explicit form of the Mellin transform (4.2). 

Since 

Ee"'^ = e"''\ /i>0, 0<v<l, (4.11) 

we have that the characteristic function ofX reads 



nv r nv 

cos — I 1 — sgn p tan — 



(4.12) 



Some manipulations as shown in D'Ovidio and Orsingher [201 1 ] prove that 



, 1 fl-r]\ 1 
EX^-i = -r — -. (4.13) 

V y V J r(i — 17} 

This permits us to conclude that 

E [N^M] "-' = e^K vr(^) nhr^) _ (4,14) 

Remark 4.3. When Xj are i.i.d Bernoulli random variables of parameter p, we have that the fractional moments of 
the compound process can be written as 

E[iV^(t)]'' = e-^^^i-P\ (4.15) 

which do not depend on r\. It follows that the mean value and the variance are 

EiV^(t) = e-^'^i-P\ VJV„(t) = e-^'(i-P^ (l - e-^'^^-P^) . (4.16) 

Note how the mean value and the variance formally coincide with those of a linear pure death process with a single 
progenitor 

Remark 4.4 (General case). Consider an infinitely divisible random variable Y in the sense of Mellin (or log-infinitely 
divisible), thus decomposable in product of i.i.d. random variables (^j. For Y we have that 



EY"-^ = [eC'^Y . (4.17) 
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W(t), 



The Mellin transform of the random product Y\ =i Cj is therefore 



E 



■Nit) 






r,-l 



Let now 9 : N ^ N such that for each r e N, 0(r) = rii=i ?j- -(f f^e random variables ^jS take integer values, then 

00 00 

E [e(JV(t))]'^-^ = XI XI m''"^Me(r) = m}Pr{JV(t) = r} (4.19) 

m=0 r=0 

00 

= 2 {e^^-^Y Pr{N[t) = r} = e-^^^+^fM^-\ 
IfQ^r) is absolutely continuous the calculations follow in the same way and arrive at the Mellin transform (4.19): 

00 fOO 

(4.20) 



E [eCJVCt))]"-^ = Y, xi-'^Pr{eiir) e dx}Pr{JV(t) = r] 

r=oJo 

oo 

= Y, {v.^'^-'y Pr{N[t) = r} = e-^P^+^^KS--^ 



r=0 



In conclusion we have that 



Nit) 



e(iv(t)) = ]~[?^-. 



(4.21) 



J=l 



4.2 Continued fractions of Cauchy random variables with Poisson distributed levels 

We consider in this section the random variables defined as 

1 



[^U^lJ ■ ■ ■ '^N(t)i — ^1 + 



x, + . 



(4.22) 



'•+'^jv(t)-i+v7r 



where Xj, j > 1, are independent Cauchy random variables with scale parameter equal to unity and location 
parameter equal to zero. We will write X ~ C(0, 1). Furthermore N[t), t > 0, is a homogeneous Poisson process 
independent of the Cauchy random variables Xj. For the convenience of the reader we note that 



[Xi]=X„ 
[X,;X2]=X, + ^, 

[X,;X2,X^]=X,+ 



X2+T-' 



(4.23) 
(4.24) 

(4.25) 



The standard Cauchy random variable has the remarkable property that X ~ 1/X, and this is the reason for 

which continued fractions can be treated when Cauchy random variables are involved (Cammarota and Orsingher 

[2010]). 

For our analysis, we need the following result. 
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Lemma 4.1. For a Cauchy random variable C(a, b), a e R, b e W^, the random variable 1/C(a, b) ~ C(a/(a^ + 
b^), b/Ca^ + b^)). In our case, a = Oand therefore 1/C(0, b) ~ C[a, 1/b). 

Our first result is stated in the next theorem. 

Theorem 4.2. The nth level fraction 

[X,;X2,...,XJ=X, + - , (4.26) 

X2 H 

""•+^„-i+3r 
has Cauchy distribution with scale parameter b„ = F^^i/F^, where F„ are the Fibonacci numbers. 
Proof. We proceed by induction. 

[X,;X2] = Xi + — ~ C(0, 2). (4.27) 



^2 



In view of Lemma 4.1, we have that 



[X,;X2,X^]=X, + ^~C(0,3/2). (4.28) 

^2 + ^ 



Furthermore, 

1 

In general we have that 



[X,;X2,X3,X^] =X,+ ,^ .^ ^, ~ C(0,5/3). (4.29) 



= Xi + C(0,F„/F„_i) = C(0, 1 + F„/F„_i) = C(0,F„+i/Fj, 
and in the last step we took into account the definition of Fibonacci numbers. D 

Remark 4.5. The Fibonacci numbers can be written in terms of the golden ration ^ = (1 + V^)/2 as 

(/,"-(l-0)" 



(4.31) 
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Therefore 

F„+i <^"+i-(l-<^)"+i ^"(^)" 

This means that [Xi;X2, ■ . . ,X^] —> C(0, (p). 

Remark 4.6. From the analysis above, we infer that [Xi;X2, ■ . . ,XjY(t)], t > 0, is a process and, for each t, possesses 
distribution equal to 

M[X,;X,,...,X.,,]).dxl/dx = i:i^^,^J±^ (4.33) 
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We give now an alternative representation to the process [Xi;X2, . . . ,Xjv(t)], t > 0. 

Theorem 4.3. The characteristic function of the random continued fraction [Xi;X2,...,Xf^(^^^] reads 

]Eg,/Jft;X„...^„„]^g-|/3|0^]-[g-l/J|^( — J pr{M[t) = n}. (4.34) 

n=0 ;"=! 

Proof. In view of Theorem 4.2, we have that 

00 

n=0 

00 ^ 

= ^e"'^'^Pr{JV(t) = n}. 



Since 



we have that 



n=0 



""' \' An (4.36) 



F„ 1 






* '-P^ S 



* J lU''- * 



"i 1 _ ^ .^ A 1 _ ^ A "0"+i) 



H§(W-^§m 



*i-(-^lS'^'-*^"' 



</- ;^^ cp 



j=0 V 



1-^ 





"J 



„i/3[Xi;X2,...^„] 



Ee 

n=0 

00 00 



71=0 j = l 



Remark 4.7. From (4.34) we can extract the following equality in distribution: 



00 00 , ,„j 

^f[e-l^l^(^) Pr{iV(t) = n}. 



-,e-\P\<l> 



The second term in (4.38) represents the effect of randomisation of the continued fraction. 
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D 



[Xi;X2,...,XN(,)]=C(0,</.)+/52cJo,f — ^j J. (4.38) 



Remark 4.8. The above analysis suggests an alternative representation of the random continued fraction as 

^ «(t)+ i 
[^i',^2>--->^Nit)]= _2_i ^J.w(t)' (4.39) 

where the Y'^ ;^(t) are independent Cauchy random variables with scale parameter equal to 1/Fjv(£)- Clearly F„ are the 

Fibonacci numbers. 

The equality (4.39) can be ascertained by writing the characteristic function as follows: 



00 Fn+l 



]gg'PZjj=l 'i,Nl 



JV(t) 



(4.40) 



n=0 j=l n=0 

which coincides with (4.35). 
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